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ABSTRACT 



We report the discovery and orbital solutions for two new OB binaries in the 
Cygnus 0B2 Association, MT311 (B2V + B3V) and MT605 (B0.5V + B2.5:V). 
We also identify the system MT429 as a probable triple system consisting of a 
tight echpsing 2.97 day B3V+B6V pair and a BOV at a projected separation of 
138 AU. We further provide the first spectroscopic orbital solutions to the eclips- 
ing, double-lined, 0-star binary MT696 (09.5V + B1:V), the double-lined, early 
B binary MT720 (BO-IV + B1-2V), and the double-lined, 0-star binary MT771 
(07V + 09V). These systems exhibit orbital periods between 1.5 days and 12.3 
days, with the majority having P<6 days. The two new binary discoveries and 
six spectroscopic solutions bring the total number of known massive binaries in 
the central region of the Cygnus 0B2 Association to 20, with all but two having 
full orbital solutions. 

Subject headings: techniques: radial velocities — binaries: general — (stars:) 
binaries: spectroscopic — (stars:) binaries: {including multiple) close — stars: 
early-type — stars: kinematics — surveys 



Introduction 



The ubiquity of close massive binaries in open clus ters i s becoming well docu mented 



(e.g.. lSana et al.ll2nn9l . iMahv et al.ll2nn9l . ISana et al.ll2nn8l . and lPe Becker et al.ll2nn6l ). How- 
ever, the number of orbital solutions for massive binaries is still relatively small. In our 
previous works docurnenting the number and statistics of massive binaries in Cygnus 0B2 
kiminki et aP 120071 boosl . l2009h . we discovered nine systems (MT059, MT145, MT252, 
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tional systems (A36 and A45; notation from IComeron et al.l |2002| ) as part of our ongoing 



radial velocity survey in one of the Galaxy's super OB clusters. The primary objective of 
the Cyg 0B2 radial velocity survey is to aid our understanding of numerous open questions 
involving massive stars, including their formation, the predicted rates of supernovae type 
Ib/c arising through binary channels, and whether there is a correlation between cluster 
density and the binary fraction therein. 

In this work, we report discovery of two new double-lined spectroscopic binaries (SB2s), 
and we provide full orbital solutions (periods, mass ratios, eccentricities and probable in- 
clinations) for these and four known binaries previously lacking full solutions. This brings 
the number of massive binaries having full or partial solutions to 20, the na ost of any single 
cluster or association. In a companion paper ( iKiminki &: KobulnickyI 120121 ). we summarize 
the current state of massive binaries in the literature and use the combined binary orbital 
parameter information and ephemerides for 113 massive stars in Cyg 0B2 to infer the intrin- 
sic orbital parameter distributions and the binary fraction of the Association. We present 
these findings as a com plement to analyses of massive sta rs in other cluste rs (such as the 
work with Sco B2 bv Kouwenhovenl boO?!. NG C 6231 bv ISana et aPboOsl . NGC 2244 by 
Mahy et al.ll2009l . and NGC 66 11 bvlSana et al.ll2009f). For a more comprehensive discussion 
of the scope of this project see iKiminki fc KobulnickyI (j2012[ ). 



Section 2 reviews the observations of this ongoing 12- year radial velocity survey and 
discusses the method we used to determine radial velocities. Section 3 discusses the method 
for computing orbital solutions and provides new spectroscopic orbital solutions for MT311, 
MT429, MT605, MT696, MT720, and MT771. Finally, Section 4 summarizes the survey 
findings to date. 



2. Observations and Radial Velocity Measurements 



Papers I-III (IKiminki et al.l 120071 . IKiminki et al.ll2008l . and IKiminki et al.l l2009l respec- 
tively) detail the observations of this survey through 2008 October (see Table 1 of Paper 
III). Since then, we have obtained eleven additional datasets with the WIRO-Longslit spec- 
trograph on the Wyoming Infrared Observatory (WIRO) 2.3 m telescope. The observation 
dates, including Heliocentric Julian Dates, for these datasets are listed in Table HJ We used 
the 1800 1 mm~^ grating in first order and obtained a typical spectral resolution of ~1.5 A 
FWHM across the chip. These observations took place over 26 nights between 2008 De- 
cember 17 and 2010 October 27 to primarily examine the Ha, He I, and He II absorption 
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lines in suspected SB2s. Exposure times for all observations varied from 300 s to 4500 s (in 
multiples of 300-900 s) depending on weather conditions and yielded a maximum signal-to- 
noise ratio (S/N) of 200:1 for the brightest stars. The spectral coverage was 5250-6750 A. 
Copper- Argon lamp exposures were taken after each star exposure to wavelength calibrate 
the spectra to an rms of 0.03 A (1.4 km s~^ at 6400 A). Spectra were Doppler corrected 
to the heliocentric frame. All new datasets were reduced using standard IRAF0 reduction 
routines as outlined in Paper I. 

We obtained radial velocities for the single-lined spectroscopic binary (SBl), MT 429 via 
the IRAF cross-correlation t ask XCSAO in the RV SAO package faurtz et al. 1992), using 
a model stellar atmosphere (ILanz fc Hubenyl |2003| . TLUSTY) of the appropriate effective 
temperature and gravity (Paper I). Wavelength regions used include AA4020-4500 A (ex- 
cluding the DIB at A4428 A) for the WIYN spectra, AA5250-6750 A (excluding all regions 
except He II A5411 A, He I AA5876, 6678 A, and Ha) for the 2008 and 2009 WIRO data. 



and A A4200-6000 A (excluding the DIB at A4428 A) for the 2010 WIRO data. Kurtz et al. 
(119921 ) calculate the uncertainties within the XCSAO task as 



3w 



:i + r) 



(1) 



a function of the r statistic (t he ratio of the correla tion peak height to the amplitude of 
antisymmetric noise) derived in iTonry fc David (119791 ). w in this equation is the FWHM of 
the correlation peak. 

In order to measure radial velocities for the SB2s, we fit multiple Gaussians to the He I 
AA5876, 6678 A and Ha lines. We obtain rough fits to the most deblended (and highest 
S/N) spectra using the SPLOT routine in IRAF. These values are then used to fit two four- 
parameter Gaussians simultaneously to those same spectra using the IDL multiple Gaussian 
fitting routine ARM_MULTGAUSS (written by Andrew R. Marble at the University of 
Arizona) , which rn akes use of the non-linear, least-squares, curve fitting IDL package MPFIT 
(IMarkwardtl l2009l ) . We then adopt the average of the Gaussian widths and depths as fixed 
parameters to fit line centers for all deblended spectra of the system. The radial velocities 
are then calculated from the best-fit Gaussian centers. Gaussian parameter uncertainties are 
obtained from the square-root of the diagonal elements of the covariance matrix multiplied 
by the measured of the fit. The overall radial velocity uncertainties are obtained by 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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summing the uncertainties on the hne centers in quadrature with the typical uncertainty in 
the wavelength calibration (1.4 km s~^ for WIRO and 2 km s~^ for WIYN). The wavelength 
solutions of the WIRO-Longslit spectra have been adjusted to a common zero point using 
the mean of the measured radial velocities of the interstellar Na I AA5889.95, 5895.92 A 
lines, which we assume are invariant for a giv en star. WIYN and Keck spectra have been 
calibrated using radial velocity standard stars ( iKiminki et al.l 120071 . 120081 ). 



3. Spectroscopic Solutions 



To estimate the periods from the radial velocity data for both the SBls and SB2s, we 
an alyze the pow e r spec trum as filtered through the the CLEAN deconvolution algorithm 
of [Roberts et al.l (119871 ). For computation of the orbital solutions w e used th e Bina ry Star 
Combined Solution Package (referred to as BSCSP in this work) of iGudehusI ( 1200 ll ) which 
makes use of a nonlinear, least-squares subroutine (NONLINFIT) to solve for orbital param- 
eters of binary (SBl and SB2) and multiple systems using spectroscopic and/or astrometric 
data. As an added check, we also computed sin gle-component spectroscopic solutions with 
the least-squares, curve fitting program SBCM (IMorbey &: Brosterhuslll974f ) and combined 
component solutions with SBTWOS, our own MATLAB-based orbital solution program uti- 
lizing Markov-Chain Monte Carlo techniques. Solutions from all programs were in good 
agreement, but BSCSP and SBTWOS yielded orbital parameter uncertainties up to three 
times larger than those from SBCM. In this work, we present the results with the more con- 
servative uncertainties f rom BSCSP. For the eclipsing, spectroscopic, triple system MT429, 
we also used P HOEBE (IPrsa fc Zwitten 120051) . an eclipsing binary modeling program based 
on the code of IWilson fc Devinneyl (jl97ll ). PHOEBE provides orbital solutions for eclipsing 
binary systems with photometric and/or spectroscopic data. 



3.1. MT311 

At V~14, the SB2 MT311 is among the faintest in our sample, and the resulting spec- 
tra have an average s i gnal-t o-noise ratio (S/N) of ~60. The primary is classified a B2V by 



Massey fc Thompson! (1199 ll ). Given the lower S/N and velocity separation of ~210 km s~^, 
we can only determine approximate spectral types for the components. An absence of He II 
and relatively strong He I places both components in the early B range. The 2001 August 24 
spectrum, obtained at WIYN, shows weak Mg II A4481 A and absent or weak silicon, nitro- 
gen, and carbon (i.e., C III A4070 A, C II A4267 A, Si IV AA4089,4116 A, Si II AA4128,4130 
A, O II AA4415,4417 A, N II A3995 A). This suggests that both components are likely still on 
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the main sequence. The temperature-sensitive ratio He I A4471 A:Mg II A4481 A indicates an 
approximate temperature class of B2-B3 for both components. Figured] shows the similarity 
between the line depths for He I A5876 A, Ha, and He I A6678 A, suggesting similar lumi- 
nosities and providing additional support for comparable temperature classes. We, therefore, 
tentatively classify the components a B2V and a B3V. In Figure [21 we show the 2001 August 
24 spectrum of MT311 obtained at WIYN (top). We provide a c omposite spectrum be 



l ow M T311, composed of HD42401 (B2V) and HD74280 (B3V) from lWalborn fc Fitzpatrick 



( 1l990l ). HD42401 and HD74280 have been shifted to the appropriate component radial ve- 



locities from Table |2] and scaled to the appropriate luminosity ratio f or a B2V and B3V. The 



MT31 1 spectrum is smoothed to approximately the resolution of the lWalborn &: Fitzpatrick 



( 1990 ) spectra. The stronger absorption present in the He I A4471 A and Mg II A4481 A 
lines in the MT311 spectrum may be a normalization effect given that the lines fall near the 
end of the chip. The remaining lines are consistent with the composite spectrum. 

Table [2] contains the ephemeris for this system, including radial velocities for the 11 
highest S/N spectra. The majority of the spectra were obtained from 2009 June 27 through 
2009 September 4 with the Wyoming Infrared Observatory (WIRO) Longslit spectrograph. 
One observation is from 2001 with the WIYn| 3.5 m telescope and Hydra spectrograph. 
With these 11 spectra, covering a span of eight years, we calculate radial velocities from 
the He I A5876 A line (He I A4471 A line for the WIYN spectrum) and obtain a period of 
5.75 days from the power spectrum for both components with no significant aliases. The best 
orbital solution, displayed in Figured yields a period of F = 5.752 ±0.005 days with a fixed 
eccentricity of zero. Allowing the eccentricity to vary results in an insignificant decrease 
in the rms of the fit and yields e = 0.22 ± 0.05. Given that the 11 data points do not 
provide a strong constraint on the eccentricity, we adopt the solution with the fixed zero 
eccentricity. In the figure, solid symbols correspond to the radial velocities of the primary, 
open symbols correspond to the radial velocities of the secondary, and the different shapes 
represent the different instruments/observatories (i.e, squares represent observations taken 
with WIRO-Longslit and triangles represent observations taken with Hydra at WIYN). We 
found it more difficult to fit the secondary's spectral features owing to the low S/N and 
heavy blending of profiles and thus only utilized 10 measurements (the measurement at 
(j) = 0.808 was excluded). With a period of F = 5.752 ± 0.005 days, an eccentricity of 
e = 0.0, and semi-amplitudes of Ki = 88 ± 5 km s^^ and K2 = 118 ± 9 km s~^, we 
calculate masses of Misin^i = 3.0 ± 0.7 M© and Mssin^^ = 2.2 ± 0.4 Mq, and semimajor 
axes of oisini = 10.0 ± 0.6 Kq and a2sim = 13.40 ± 0.01 R0. This yields a mass ratio of 



^The WIYN Observatory is a joint facility of the University of Wisconsin, Indiana University, Yale 
University, and the National Optical Astronomy Observatory. 
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q = 0.7 5 =1= 0.07 and is consistent with the theoretical mass ratio of g = 0.76 for a B2V+B3V 
system (jPriUing fc Landoltll200Cll ) . The calculated lower limit masses indicate that the system 



needs to be inclined at i ~42° to be consistent with the theoretical masses of 10 and 



7.6 Mq for a B2V star and B3V star flDrilling fc Landoltl 120001 ). 



The computed heliocentric systemic velocity for MT311 is 7 = 1±4 km s~^. This differs 
by nearly 11 km s~^from the average systemic velocity of —10.3 km s~^ for other members of 
the Association (Pa per I), calling into question its membership. If we adopt an apparent V 
magnit ude of 13.87 (IMassev fc Thompson! 1 199 111 for the system, the B-V color of 1.39 for the 
system ( iMassey &: Thompsorull99l(). an absolute V m agnitudes of -2.45 and -1.94 for a B2V 
and B3V respectively f lHumphreys fc McElroylll984l ). an intrinsic color [B — V)^ = —0.21 
and a reddening law with i?y=3.0, the resulting spectrophotometric distance to MT311 is 
~2.6 kpc (dist ance modulus of 12.07). The accepted distance to Cyg 0B2 is near 1.7 kpc 
( lHansonll2003l ) but with a relatively large uncertainty owing to the uncertainty on absolute 
mag nitudes for massive sta rs and the observed variability in the total-to-selective extinction, 
Rv (jPatriarchi et al.l 120031 ). Although our computed spectrophotometric distance is 50% 
larger than the accepted distance to Cyg 0B2, the uncertainties are sufficiently large that 
we cannot draw a definitive conclusion regarding the membership of MT311. 

Table [2] provides the ephemerides for MT311, MT605, MT696, MT720, and MT771, 
listing the date (HJD— 2,400,000), phase (0), measured radial velocities (Ki and K2), and 
observed minus calculated velocities (Oi — Ci and O2 — C2). Error estimates are shown 
in parentheses. The complete list of orbital elements and physical parameters for MT311, 
MT429, MT605, MT696, MT720, and MT771 appears in Table l Listed within the table 
are the period in days (F), eccentricity of the orbit (e), longitude of periastron in degrees 
(oj), systemic radial velocity (7), epoch of periastron (Tq), primary and secondary semi- 
amplitudes {Ki & K2), adopted or calculated minimum primary and secondary masses in 
solar masses (Mi & M2), primary and secondary mass functions in solar masses (/(m)i & 
f [171)2), spectral classifications from this survey (S.Ci & S.C.2), the minimum primary and 
secondary semi- major axes in solar radii (aisin i & a2sin i), and finally, the rms of the best 
fitting orbital solution [rmsi Sz rms2)- Additionally, the inclination, primary and secondary 
effective temperatures, and component stellar radii are included for MT429. 



3.2. MT429 



iMassev fc Thompson! f ll99l[ ) classify MT429 (=V2186 Cyg) as a BOV. IPigulski fc Kolaczkowski 
( 119981 ) showed that this is an eclipsing system of the Algol type, estimating a period of 2.9788 
days based on an incomplete light curve. Our spectra of MT429 initially suggested it to be 
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an SBl with low-amplitude radial velocity variations at a very similar period. 

Figure H] shows a slightly smoothed, high quality spectrum of MT429, obtained 2001 
September 8, with 09V (HD46202), BOV (H D36512), and BIV (HD144470) spectra from the 



digital atlas of lWalborn fc Fitzpatrick for comparison. With the exception of the lower 



S/N 2000 September 18 Keck spectrum, our spectra of MT429 do not cover the important 
Si III A4552 A line, used with the Si IV A4089 A line and the Si II AA4128-4130 A lines as 
the primary temperature class indicator in early B stars. However, the WIYN spectrum in 
the figure shows a moderately weak Si IV A4089 A line and a moderately strong C III + O II 
blend at AA4070 A. These lines qu ickly diminish with decreasing temperature in early B 



stars (IWalborn fc Fitzpatricklll990[ ). There is also no evidence for He II A4200. This agrees 
with a temperature class of BO-Bl. The primary luminosity indicator (Si IV AA4089,4116 A 
to the nearby He I lines) also agrees with a V luminosity class. All of our spectra support 
the classification of MT42 9 as B OV, within a fraction of a temperature class, consistent with 



the iMassey fc Thompson! ( 1l99ll ) classification. 



We obtained 16 spectra on MT429 from 2000 September 18 through 2010 October 27, 
including one at Keck, five at WIRO, and 10 at WIYN. The WIYN spectra were cross- 
correlated using a wavelength range of AA4020-4500 A (excluding the DIB at A4428 A). 
The WIRO 2010 spectra were cross-correlated using a wavelength range of AA4200-6000 A 
(excluding the DIB at A4428 A and the interstellar (IS) Na I lines at AA5890,5895 A). The 
2008-2009 WIRO spectra were cross-correlated using a wavelength range of AA5400-6700 A 
(excluding all but the He I, He II, and hydrogen lines owing to the large number of IS lines 
and DIBs in this spectral range). The Keck spectrum was cross-correlated using a wavelength 
range of AA4050-5200 A (excluding the DIB at A4428 A) 

The strongest peak in the power spectrum corresponds to a pe riod of P = 2.9786 days 
with n o sig nificant aliases, consiste nt with the photometric periods of lPigulski fc Kolaczkowski 



( 119981 ) and [Henderson et al.l (1201 if ). Allowing all orbital parameters to vary, BSCSP provides 
an orbital fit with P = 2.9785, e = 0.24, K = 17.7 km s"\ and rms of 3.0 km s'^ We 
excluded one outlier from the fit, the 2008 June 15 WIYN measurement {Vr = —48 km s^^). 
Fixing the eccentricity at zero provides a comparable fit to the data with a similar rms of 
(3.4 km s~^), period of P = 2.9786±0.0001 days and systemic velocity of 7 = — 14±1 km s^^. 
The resulting mass function of /(m) = 0.0012 ± 0.0003 would suggest that either the sec- 
ondary is a low-mass star (Mi ~ 0.8 M0) or that the inclination is very small. A small 
inclination is inconsistent with the observed eclipses. A large inclination coupled with a 
low-mass companion could conceivably explain the observed radial velocity curve, but such 
a scenario entails that only primary eclipses would be observed given the vast difference 
in temperature between a main sequence or evolved ~1 MQsecondary and a BOV primary. 
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However, existing photometric data clearly show both a primary and a secondary eclipse. 



We combined the photometric data of iPigulski fc Kolaczkowskil (119981 ) with that of 
Henderson et al.l (120111 ) who also obtained I-band photometry, finding a nearly identical 
period of 2.9788 days. We adopted an apparent I magnitude of 11.08 based on the mean 
maxim um magnitude from [Henderson et al.l (120111 ). There are 3,7 10 days between the end 
of the iPigulski fc Kolaczkowskil (119981) data and th e start of the [Henderson et al.l (120111 ) 
data. This necessitates shifting the [Henderson et al. ( 2011 ) data by phase = 0.225. This 
is approximately equal to a 0.0002 day difference in the calculated period but is consis- 
tent within the calculated period uncertaintyB Figure |5] shows the combined, adjusted 
photometric data folded at the 2.9786 day spectroscopic period. Open circles den o te the 
Pigulski &: Kolaczkowskil (1l998l ) data, and the crosses denote the iHenderson et al.l (l201l[ ) 
data. The data show a primary to secondary eclipse depth of slightly less than 2:1, indi- 
cating components with similar temperatures. Adopting the period, eccentricity, epoch of 
periastron, and longitude of periastron from the spectroscopic solution, the data were best 
fit with an inclination of i = 74.3°, mass ratio of q = 0.59, separation of asini = 22.5 R©, 
primary radius of Ri = 5.5 R©, secondary radius of R2 = 3.9 R©, primary temperature of 
Ti = 20, 900 K, and secondary temperatur e of T2 = 15, 200 K . These results are consistent 
with a B2V primary and B5V secondary iDrilling &: Landoltl (I2OOOI ). However, the semi- 
amplitude of 16 km s~^ measured from the spectroscopic solution does not agree with the 
predicted 137 km s"^ primary radial velocity semiamplitude from PHOEBE. 

The seemingly inconsistent conclusions from the light curve analysis (a B2V + B5V 
eclipsing binary) and the spectroscopic velocity curve analysis (a BOV with a velocity am- 
plitude of 16 km s~^ implying a very low-mass com panion) posed a conundrum until the 
the discovery that MT429 is a close visual double. ICabellero-Nieves et al.l ( 201ll ) present 
adaptive-optics K-band images and HST fine guidance sensor measurements showing that 
the MT429 system is a visual double with a separation of 0'.'08 (138 AU projected separation 
at a distance of 1.7 kpc). The fainter component, MT429b, is 1.05 mag fainter at V band 
than MT429a, consistent with an early-to-mid B star and a BOV at a common distance and 
reddening. A consistent picture of MT429 emerges if it is actually a triple system consist- 
ing of a BOV star (MT429a) separated by at least 138 AU from a tight B2V-I-B5V binary 
(MT429b+MT429c) zn a 2.9786 day orbit. Statistically, the small angular separation of 0'.'08 
between the two dominant components in MT429 is improbable unless there is a physical 
association ( ICabellero-Nieves et al.l 1201 ll ). With this new information, we adjusted the pho- 
tometric data to remove the additional flux contributed by the BOV, estimated to be ~72% of 



•^The radial velocity data also required a phase shift amounting to 0.09 after 1,137 days, which is equivalent 
to a period difference of ~0.0002 days and consistent within our calculated orbital period uncertainty. 
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the total flux by lCabellero-Nieves et al.l ( 1201 ll ). Because a larger percentage of the total light 



is being eclipsed with the BOV flux subtracted, the eclipse depths appear deeper compared 
to the hght curve in Figure O The upper panel in Figure M shows the adjusted, combined 
photometric data with the new PHOEBE fit. Again, adopting the period, eccentricity, epoch 
of periastron, and longitude of periastron from the spectroscopic solution, the best fit to the 
data results in an inclination ofi = 89°, mass ratio of g = 0.68, separation of asini = 20.4 Rq, 
primary radius of Ri = 4.8 R©, secondary radius of i?2 = 3.6 Rq, primary temperature of 
Ti = 19, 100 K, and secondary tempe rature of T2 = 13, 700 K. This is most consistent with 



a B3V primary and B6V secondary (jPrilling &: Landoltl I2OOOI ) . The new result predicts a 



radial velocity primary semiamplitude of ~140 km s~^ (205 km for the secondary), still 
in sharp disagreement with the measured semiamplitude of 16 km s~^. The lower panel of 
Figure M shows the theoretical velocity curves of the B3V (solid curve) and B6V (dashed 
curve) stars prescribed by the above scenario. We believe the answer to this discrepancy lies 
in the fact that the B3V is roughly 5.5 x fainter than the BOV and the B6V is significantly 
fainter than the B3V in the visual band, and the large velocity amplitude of the tight bi- 
nary is masked (spectral dilution) by the much brighter (apparently) constant-velocity BOV 
component MT429a. Simulations wherein we combine spectra with appropriate flux and 
signal-to-noise ratios support the conclusion that a 140 km s~^ velocity amplitude signal 
from MT429b is diluted to the level of ~16 km s~^ in a composite B0V-I-B3V spectrum. 
The symbols in the lower panel (filled circles for WIRO, triangles for WIYN, a diamond for 
Keck) show the apparent velocity curve measured in the composite spectrum of the system. 
The dotted line shows the systemic velocity of the binary (7 = —14 ± 1). While the period 
of the folded velocity data agrees well with the photometric data, the small implied velocity 
amplitude of 16 km s~^ is not representative of the velocity of any component within the 
system. MT429, then, is the first massive triple system discovered in the CygOB2 radial 
velocity survey. It joins Schulte 5 and Schulte 8 as a growing group of N>2 multiple-star 
systems in the Cyg 0B2 Association. 

We adopt from BSCSP the values and uncertainties pertaining to the orbital period, 
eccentricity, epoch of periastron, and angle of periastron for the B3V+B6V pair and list 
these in Table |31 We adopt the values for component radii, masses, effective temperatures, 
theoretical semiamplitude semimajor axes, and mass ratio from PHOEBE. These pertain 
to the tight B3V+B6V system. The listed systemic velocity of the system pertains to the 
composite system. The RV ephemeris for MT429 is listed in Table H] and includes the date 
{HJD — 2,400,000), phase (0), measured radial velocity (K ), cross-correlation la error 
{la err). We reiterate here that, while these velocities show excellent agreement with the 
period and phase angle derived from the light curve, the amplitudes are not physically mean- 
ingful, as they reflect the radial motion of the B3V component as diluted by the (apparently 
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constant velocity) spectrum of the brighter BOV component. 



3.3. MT605 



Massey fc ThompsonI (1l99lh classify MT605 a B0.5V. In Paper I, we typed the presumed 



single system a BIV based on visual comparisons with spectra from the lWalborn fc Fitzpatrick 



( 119901 ) digital atlas. We now know that this is an SB2. Examination of the highest quality 
spectra for this system shows that neither component's spectrum displays any He II, indi- 
cating that the components are likely later than BO. However, very weak Mg II A4481 A 
absorption relative to He I A4471 A indicates the components are earlier than B2. Further- 
more, because Si III A4552 A and Si IV A4089 A are also weak and C III+O II AA4070, 
4650 A absorption is strong, we reclassify the primary a B0.5V (in agreement with the origi- 



nal assessment of iMassey fc ThompsonI Il99ll ) with acknowledgment that the small velocity 



separation (~160 km s~^ for He I) at a resolution of R~4500 makes the luminosity classifi- 
cation difficult and still somewhat uncertain. With a mass ratio near unity and a luminosity 
ratio also near unity (determined from examination of all lines while in or near quadra- 
ture), we also tentatively classify the secondary a B0.5-B1 with an indeterminate luminosity 
class. Figure [7] displays the 2001 September 9 spectrum of MT6 05 taken near quadra- 



ture ( top) and obtained at WIYN and HD36960 (B0.5V) from the IWalborn fc Fitzpatrick 



( ll990l ) digital atlas for comparison . The WIYN spectrum has been boxcar smoothed to 



approximately the resolution of the IWalborn fc Fitzpatrickl ( 1l990l ) spectra. The comparison 



spectrum is in good agreement with the WIYN spectrum of MT605, with the exception of 
the C III+O II AA4070,4650 A blend and Si IV A4089 A absorptions, which appear shghtly 
stronger in HD36960. 

We observed this system 37 times between 1999 and 2009, obtaining 10 spectra (all 
with WIRO-Longslit) with sufficiently high enough S/N and large enough component profile 
velocity separations (~150 km s~^ at max) to measure radial velocities. Because of the small 
number of suitable spectra available for measuring radial velocities, determining an estimate 
of the period from the power spectrum is difficult. We, therefore, fit single Gaussians to the 
He I A5876 A and Ha lines and measured the profiles' Gaussian FWHM in 18 of the WIRO 
spectra (17 for He I A5876 A). From these measurements, we observe a strong peak in the 
power spectrum at ~6 days for both He I and Ha. This represents half the period of the 
binary and tells us the location of the primary and secondary spectral features when the 
profiles are blended. Using this information, we performed the same technique as with the 
other SB2s to calculate radial velocities from the He I A5876 A line for the 10 spectra. The 
power spectrum revealed that the likely period from the 10 spectra is 12.26 days for both the 
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primary and secondary star spectral features. The best orbital solution provided by BSCSP 
and SBTWOS yields a period of P = 12.27 ± 0.02 days, an eccentricity of e = 0.24 ± 0.07, 
and semi-amplitudes of Ki = 44 ± 3 km s~^ and = 77 ± 5 km s~^. This solution is 
shown with the folded radial velocity curve in Figure [Si The solution also results in masses 
of Misin^i = 1.3 ± 0.3 and Mssin^z = 0.7 ± 0.1, indicating a mass ratio of 13' = 0.57 ± 0.05. 
This computed mass ratio suggests that the secondary may be 2-3 subclasses later than the 
primary, and a main-sequence spectral type of B2.5:V would be in better agreement with 
the observed mass ratio than our original determination of a Bl from analysis of the line 
profiles alone. In addition, a relatively low inclination of z ~ 26° would be required for 
the component masses to be consistent with the theoretical masses of a B0.5V and B2.5V 
dPriUing fc Landoltl[200ol ). 

A time series of the 10 spectra used in the analysis, in order of phase, is provided for 
illustrative purposes in Figure [H Panel 1 shows the line profiles of He I A5876 A panel 2 
shows the line profiles of Ha, and panel 3 shows the line profiles of He I A6678 A. 



3.4. MT696 



Rios fc DeGioia-EastwoodI ( l2004l ) reveal the binary status of this interacting, eclipsing 
system composed of an 09.5V and an early B star. The authors find a period of P = 
1.46 days and present spectra also showing the SB2 nature of the system. We present the 
first spectroscopic solution for this binary using 17 spectra obtained with WIRO-Longslit 
from 2008 August 20 through 2008 September 19. Owing to the large profile widths of the 
He I A5876 A lines and a maximum line- velocity separation of ~550 km s~^, the positively 
shifted He I A5876 A profiles and the interstellar Na I A5890 A line are sometimes partially 
blended. The broad spectral features, present in both components and possibly indicative 
of large rotational velocities, is shown in Figure [TOl The figure presents a time progression 
of the He I A58 76 A, Ha, and A6678 A lines (left , middle, and right panels respectively) in 



order of phase. iRios fc DeGioia-EastwoodI (|2004| ) report this system as W UMa type, but 



the spectra show no evidence of emission lines which are common in contact binaries and 
systems with accretion disks. Given that this system has such a short period, the possibility 
that this is (or has been) a contact binary cannot be ruled out, but the early spectral types of 
the components, the current absence of emission features, and the longer than 1 day period 
suggest that this may in fact be a Beta Lyr type system instead. 

We computed orbital solutions based on both the He I A5876 A velocitie s and Ha veloci- 



ties. T he He I and Ha spectroscopic periods for both components agreed with lRios fc DeGioia-Eastwood 



( I2OO4J ) and showed strong signals in the power spectra at ~1.46 days. The orbital solu- 
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tion based on the Ha velocities, illustrated by the solid line in Figure [TTl yields a pe- 
riod of P = 1.4692 ± 0.0005 days, eccentricity of e = 0.01 ± 0.02, semi-amplitudes of 
Ki = 261 ± 4 km s'^ and K2 = 277 ± 6 km s-\ and mass ratio of g = 0.94 ± 0.02. The Ha 
radial velocities are represented by the open (secondary) and filled (primary) circles, and the 
He I A5876 A velocities are represented by open and filled triangles. The orbital solution 
based on the He I A5876 A velocities yields a period of F = 1.467 ± 0.004 days, eccentricity 
of e = 0.05 ± 0.02, semi-amplitudes of = 261 ± 6 km s'^ and K2 = 327 ± 10 km s-^ 
and mass ratio of q = 0.80 ± 0.03. The period, eccentricity, epoch of periastron, longitude 
of periastron, and systemic velocities for both solutions showed good agreement. We find a 
systemic velocity of 7 = — 10 ± 3, consistent with other Cyg 0B2 members. The secondary 
semi-amplitudes, however, differ by ~50 km s"^. The more positive semi-amplitude of the 
Ha line may be a result of some interaction between the two stars as Ha is more easily dis- 
torted in such scenarios because of its formation higher in the stellar atmosphere or possibly 
in an extended photosphere/envelope. However, given that the Ha line profiles are high 
S/N, we adopted the period, eccentricity, systemic velocity, angle of periastron, and epoch of 
periastron from this solution and fixed them in the determination of a revised He I A5876 A 
solution. This is represented by the dotted line in Figure [TTl The ephemerides for both the 
Ha solution and the combined Ha+He I A5876 A solution are provided in Table El The 
orbital elements to the combined Ha-|-He I A5876 A solution are provided in Table [31 

Adopting the revised He I orbital solution, we obtain a mass ratio of 0.85 ± 0.03. This, 
coupled with the nearly equal line depths for all features, suggests a secondary spectral 
type similar to the primary's. Without the photometric data, we can only speculate on 
the incli nation of MT696. However, assuming a theoretical mass of 16±0.5 Mq for the 



primary ( iMartins et al.ll2005l ). the calculated masses (Misin i = 15.1±0.7 Mq and M2sm i = 
12.8 ± 0.5 Mq) suggest that the system has an inclination of near 80° and a secondary with 
a spectral type of BIV. 



3.5. MT720 

In Paper II, we find that the SB2 MT720 is most likely composed of two early B stars 
based on the absence of Mg II A4481 A and strength of the He I lines. Based on new 
WIRO-Longslit observations of this system, we agree with our previous assessment. The 
new spectra show relatively strong He I absorption and no evidence of He II. The paucity 
of stellar lines in this spectral regime and the low S/N of our spectra still make accurate 
spectral typing difficult. Figure [121 shows the dominant stellar lines in the 5250-6700 A regime 
(He I A5876 A, Ha, and He I A6678 A) near quadrature. The spectrum is a composite of 
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seven spectra between (j) = 0.5 and (j) = 0.7. The line depths and widths are similar and 
suggest comparable spectral types for the components. In order to draw further conclusions 
regarding the component spectral types, we examined the following orbital solution and 
component mass estimates. 

We measured radial velocities for 32 of the 35 WIRO-Longslit spectra obtained on 
MT720 using the He I A5876 A profiles. These 32 radial velocities provide a clear period of 
4.36 days from the power spectra of both components. The power spectra also contained 
no significant aliases. The 4.36 day period agrees with our original visual assessment of just 
under 5 days in Paper II. The best combined orbital solution obtained from BSCSP produces 
a period of P = 4.3622 ± 0.0003 days and an eccentricity of e = 0.35 ± 0.02. Figure [13] 
displays this solution. We note that this solution presents an unusually large eccentricity 
for such a short-period system. We computed additional solutions with an eccentricity fixed 
at lower values, but these resulted in a much larger rms (>>21.6 km s^^ for the primary 
and >>29.5 km s~^ for the secondary). Given semi- amplitudes oi Ki = 173 ± 5 km s"^ 
and K2 = 242 ± 7 km s~^, the calculated mass ratio is g = 0.71 ± 0.06. Additionally, the 
calculated minimum masses, Misin^i = 15.5 ±0.8 Mq and Masin^z = 11.1 ±0.6 Mq, indicate 
that the likely spectral types of the compon ents are BO-IV and B1- 2V for the primary 



and secondary respectively (interpolated from iDrilling fc Landoltll2000[ ). There is no record 



of photometric variability for this system, but if the primary has a BOV or later spectral 
type, the minimum masses of the components preclude inclinations lower than i ~ 70°. The 
remaining orbital parameter values are listed in Table [31 



3.6. MT771 

The He I lines of the SB2, late 0-t ype binary MT771 show a strong Struve-Sahade effect 



(for a good summary of this effect, see lLinder et al.ll2007| ). This effect causes asymmetries or 



varying line depth/width and, in the case of MT771, is present in both components' spectra. 
This makes component identification tricky for such a short period (< 3 days) system, even 
with a maximum velocity separation of ~300 km s^^. The much broader Ha line remains 
heavily blended even at the largest velocity separation. We used this to our advantage, 
following the same procedure as with MT605, fitting a single Gaussian profile to the Ha line 
and measuring the Gaussian FWHM. The period produced by the FWHM measurements is 
equal to half the binary period. The FWHM measurements yielded a strong signal in the 
power spectrum corresponding to a period of 1.43 days with no significant aliases. Even 
with the varying asymmetries in the component He I profiles, we still chose to use the 
average measured Gaussian width and depth for each component's He I A5876 A line and 
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fix these values when measuring the radial velocities (assuming a period of ~2.86 days). 
Of the 34 epochs included for MT771 in Table [2l we omitted two observations for the 
secondary {HJD = 2, 454, 700.85 & 2, 455, 041.67) when the lines were at their most blended 
state. As was expected, the power spectra for both components showed a clear, strong 
signal at 2.86 days. The computed orbital solution for this system produces a period of 
P = 2.8635 ± 0.0002 days, an eccentricity of e = 0.05 ± 0.03, and semi-amplitudes of 
Ki = 139 ± 5 km s~^ and K2 = 164 ± 6 km s~^. The solution also produces lower limit 
component masses of Misin^i = 4.5 ±0.3 M© and Masin^z = 3.8 ±0.2 Mq and suggests a low 
inclination of i 5, 37° for this systeni (the approximate inclination required to yield the mass 



of a late O star; iMartins et al.ll2005[ ). This solution and the radial velocity curve, produced 
from the 34 WIRO-Longslit spectra obtained between 2007 August 28 and 2009 August 3, 
is shown in Figure HH The radial velocities measured at (j) =0.8-1.0 appear systematically 
more positive, and this likely stems from the strong asymmetries at these epochs. This 
is more easily seen in Figure [151 which shows a time series of the He I A5876 A line in 
order of phase. Panel 1 contains the line profiles for =0.07-0.45 and panel 2 contains 
the line profiles for (p =0.45-1.00. The He I A5876 A line of the primary appears much 
stronger than the secondary at phases of (j) =0.856-0.906, and the secondary also appears 
asymmetric. The ratio of He line depth also appears to change over a small range of phase 
(e.g., (j) =0.355-0.387), though less dramatically. 

The strong Struve-Sahade effect indicates that there are likely significant irradiation 
effects in this system, possibly stemming from a shock region between the components' 
stellar winds. If so, the shock region is nearer the secondary as the asymmetry is stronger in 
the secondary's spectral lines. The asymmetry may be the result of an unresolved emission 
feature between the component lines and originate from the shock region. The seemingly 
irregular strength of the emission would indicate varying wind velocities and densities of one 
or both stars. 

The Struve-Sahade effect also makes spectral typing difficult given that the He I and 
He II lines are crucial to temperature classification in late O stars (e.g.. He II A4541 A:He I A4471 A, 
He II A4541 A:Si III A4552 A, and He I A4471 A:Mg II A4481 A). Our initial classifications 
of 07V for the primary and 09V for the secondary, however, appear to still be close given 
the new computed mass ratio of q = 0.85 ± 0.04, the moderately strong He II at AA4200, 
4542, and 5411 A, and strong He I. 
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4. Summary 

We present orbital solutions for two new SB2s, MT311 and MT605. We also provide 
the first spectroscopic solutions to the eclipsing system, MT696, the early B, eclipsing triple 
SBl, MT429, the early B SB2, MT720, and the late-0 SB2, MT771. 

Table |5] lists the 20 known OB binaries in Cyg 0B2. Column 1 gives the system 
designations. Column 2 shows the binary types (i.e., eclipsing, SBl, SB2, etc). Column 3 
provides the spectral types of the components. The spectral types in square brackets have 
been determined indirectly from the data (i.e., based on the computed mass ratio rather 
than the component's spectrum). Columns 4-6 give the orbital periods, eccentricities, and 
mass-ratios respectively. Finally, column 7 lists the corresponding surveys responsible for 
the orbital information. Figure [16] shows the locations of these systems relative to the other 
OB stars in this survey. There are no apparent clusterings or groups of massive stars, nor are 
there tendencies of the massive binaries to reside in the center of the Association versus the 
outskirts. Of the 20, we know reliably the orbital periods of all 20, the orbital eccentricities 
of 18, and the mass-ratios of all 20 (the mass ratios for the SBls are actually range estimates 
based on the range of probable inclinations). There are six SBls, 14 SB2s, five eclipsing 
binaries of the Algol type, two eclipsing binaries of the j3 Lyr type, and one eclipsing binary 
of the W UMa type. At least five binaries are bo th spectroscopic and eclipsing systems. In 



a companion paper (IKiminki fc Kobulnickyll2012l ). we use the orbital parameter information 



and ephemerides associated with the short period systems (P < 26 days) listed in Table O 
along with the ephemerides of 94 additional OB stars in Cyg 0B2, to compute the intrinsic 
binary orbital parameter distributions (i.e., period, mass ratio, and eccentricity) and intrinsic 
binary fraction of the cluster. 
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Table 1. Dates of WIRO Observing Runs 





Date Coverage 


Date 


(HJD) 


2008 Dec 17 


2,454,818 


2009 Jan 3 


2,454,835 


2009 Jan 16 


2,454,848 


2009 Jun 27-30 


2,455,010-2,455,013 


2009 Jul 5 


2,455,018 


2009 Jul 27,28 


2,455,040-2,455,041 


2009 Jul 30-Aug 3 


2,455,043-2,455,047 


2009 Aug 26-29 


2,455,070-2,455,073 


2009 Aug 31-Sep 4 


2,455,075-2,455,079 


2010 Aug 21 


2,455,430 


2010 Oct 27 


2,455,497 
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Table 2. Ephemerides for MT311, MT696, MT720, & MT771 











Vr2 


O2-C2 


Date (HJD-2,400,000) 




(km s ^ ) (km s "'") 


(km s^"'^) 


(km s~^) 


MT311 


52,146.68 


0.569 


-98.0 (10.5) 


-13.2 


104.5 (22.4) 


-11.9 


55,010.64 


0.457 


-87.1 (6.2) 


-10.4 


101.1 (11.0) 


-4.4 


55,011.78 


0.655 


-75.7 (6.8) 


-12.8 


83.8 (12.0) 


-3.3 


55,012.72 


0.818 


27.2 (6.3) 


7.6 


-62.3 (11.0) 


-39.2 


55,013.82 


0.010 


77.8 (4.9) 


-10.8 


-104.1 (8.5) 


11.1 


55,070.66 


0.891 


67.5 (6.8) 


11.4 


-123.3 (11.9) 


-51.5 


55,071.72 


0.075 


72.7 (4.7) 


-14.2 


-91.6 (8.4) 


21.2 


55,076.81 


0.961 


81.6 (6.9) 


1.3 


-91.6 (12.2) 


12.5 


55,077.74 


0.122 


86.5 (5.1) 


9.9 


-104.1 (9.0) 


-4.9 


55,078.73 


0.293 


20.8 (5.1) 


24.3 






55,079.76 


0.472 


-81.6 (5.1) 


-1.4 


132.3 (9.0) 


21.9 


MT605 


54,643.87 


0.241 


-52.8 (4.9) 


-3.7 


65.8 (6.8) 


11.8 


54,674.94 


0.776 


22.2 (5.1) 


-0.8 


-64.7 (7.2) 


4.4 


55,043.65 


0.837 


26.6 (5.9) 


-6.3 


-75.9 (8.2) 


13.3 


55,044.75 


0.926 


42.6 (7.2) 


5.4 


-107.8 (10.0) 


-11.1 


55,045.68 


0.003 


7.5 (6.3) 


-9.7 


-58.0 (8.8) 


4.0 


55,072.64 


0.204 


-47.7 (7.5) 


-1.4 


38.4 (10.4) 


-10.6 


55,073.69 


0.287 


-36.3 (7.4) 


13.3 


49.8 (10.4) 


-4.9 


55,075.62 


0.444 


-48.0 (6.1) 


-10.6 


37.8 (8.4) 


4.3 


55,078.63 


0.690 


7.1 (8.7) 


3.5 


-31.0 (12.2) 


6.9 


55,079.68 


0.776 


27.4 (7.7) 


6.0 


-77.4 (10.7) 


-8.3 




MT696 (He I A5876 A + 


Ha) 






54,699.74 


0.080 


-247.5 (5.2) 


9.5 


355.0 (7.2) 


73.7 


54,701.91 


0.556 


298.8 (4.2) 


50.3 


-305.2 (11.7) 


9.0 


54,724.62 


0.012 


-266.9 (3.7) 


1.0 


321.6 (2.8) 


27.4 


54,724.73 


0.087 


-262.1 (2.9) 


-9.4 


270.7 (2.4) 


-5.7 


54,724.82 


0.153 


-170.9 (4.9) 


25.3 


243.4 (7.9) 


33.6 


54,725.61 


0.687 


147.9 (2.7) 


15.6 


-186.4 (5.6) 


-9.1 


54,725.74 


0.773 


2.8 (12.5) 


4.4 


-11.9 (4.0) 


7.6 


54,725.85 


0.854 


-117.4 (4.5) 


10.8 


153.2 (11.1) 


23.6 


54,726.63 


0.381 


107.2 (3.3) 


-39.0 


-242.6 (8.5) 


-49.1 


54,726.74 


0.459 


217.9 (4.3) 


-9.1 


-308.9 (15.6) 


-20.1 


54,726.84 


0.526 


243.2 (1.9) 


-9.0 


-322.4 (8.0) 


-3.8 


54,729.63 


0.423 


165.4 (3.8) 


-30.5 


-283.2 (11.0) 


-31.0 


54,729.72 


0.488 


241.8 (3.0) 


-1.7 


-324.9 (11.2) 


-16.6 


54,729.81 


0.548 


264.5 (4.5) 


14.3 


-312.1 (9.8) 


4.2 
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Table 2 — Continued 







Ki 


Ol 




Vr2 


O2 - C2 


Date (HJD-2,400,000) 


<t> 


(km s-i) 


(km 


s-i) 


(km s-i) 


(km s-i) 


MT(i<)(i (lid ) 


54,699.74 


0.079 


-242.9 (4.7) 




14.4 


282.5 (12.3) 


29.4 


54,700.77 


0.781 


9.3 (6.0) 




22.6 


10.1 (10.4) 


16.3 


54,700.90 


0.869 


-183.5 (6.6) 




-33.7 


127.6 (7.8) 


-11.2 


54,701.91 


0.555 


254.1 (4.5) 




5.4 


-280.3 (6.6) 


4.2 


54,722.84 


0.800 


19.2 (19.2) 




63.2 


-32.3 (135.5) 


-58.7 


54,724.62 


0.011 


-269.9 (5.0) 




-2.1 


262.2 (7.3) 


-2.0 


54,724.73 


0.087 


-245.6 (4.8) 




7.6 


271.5 (11.0) 


22.7 


54,724.82 


0.152 


-180.7 (5.5) 




16.5 


197.7 (10.4) 


8.5 


54,725.61 


0.686 


160.3 (12.5) 




26.8 


-123.6 (12.9) 


38.5 


54,725.74 


0.773 


7.0 (13.0) 




7.3 


25.8 (54.8) 


45.8 


54,725.85 


0.853 


-145.2 (4.8) 




-18.2 


109.3 (2.1) 


-5.4 


54,726.63 


0.381 


156.7 (8.0) 




11.6 


-179.2 (11.4) 


-4.8 


54,726.74 


0.458 


222.5 (4.4) 




-3.9 


-279.8 (10.3) 


-19.0 


54,726.84 


0.526 


256.7 (8.7) 




4.5 


-291.3 (7.7) 


-3.0 


54,729.63 


0.423 


196.0 (5.8) 




1.0 


-252.8 (7.8) 


-25.3 


54,729.72 


0.487 


246.0 (9.5) 




2.8 


-280.1 (8.8) 


-1.4 


54,729.81 


0.548 


262.1 (9.5) 




11.7 


-287.2 (8.8) 


-0.8 


MT720 


54,341.83 


0.072 


-93.4 (10.5) 




0.1 


135.0 (14.8) 


25.9 


54,343.81 


0.524 


149.2 (15.9) 




6.9 


-228.0 (22.4) 


-7.0 


54,345.82 


0.986 


-66.8 (14.4) 




-14.0 


78.2 (20.3) 


26.2 


54,346.87 


0.226 


-166.8 (15.4) 




-19.3 


208.6 (21.7) 


23.8 


54,348.76 


0.660 


163.7 (13.5) 




3.0 


-216.6 (19.0) 


30.4 


54,399.56 


0.306 


-171.8 (7.1) 




-17.8 


180.5 (9.9) 


-13.3 


54,399.71 


0.340 


-153.4 (7.2) 




-7.8 


179.2 (10.1) 


-2.9 


54,401.55 


0.760 


95.8 (7.3) 




9.4 


-130.8 (10.3) 


12.2 


54,401.71 


0.798 


83.4 (8.3) 




23.5 


-96.8 (11.7) 


8.9 


54,402.62 


0.006 


-79.0 (5.8) 




-16.2 


76.8 (8.2) 


10.8 


54,402.75 


0.036 


-76.5 (7.0) 




0.6 


88.8 (9.8) 


2.7 


54,405.60 


0.691 


138.9 (9.7) 




0.5 


-195.5 (13.6) 


20.2 


54,406.59 


0.916 


-49.4 (7.8) 




-34.7 


33.9 (11.2) 


35.3 


54,408.58 


0.374 


-141.9 (9.5) 




-15.5 


128.2 (13.3) 


-26.9 


54,409.57 


0.600 


183.0 (9.9) 




-7.0 


-256.7 (13.9) 


31.2 


54,409.75 


0.641 


139.6 (9.7) 




-32.9 


-222.8 (13.7) 


40.7 


54,410.57 


0.830 


39.2 (7.2) 




0.9 


-114.3 (10.1) 


-39.0 


54,410.75 


0.871 


49.2 (11.5) 




37.3 


-62.2 (16.1) 


-23.6 


54,643.88 


0.313 


-141.7 (5.6) 




11.3 


216.9 (8.0) 


24.4 


54,696.75 


0.434 


-70.8 (6.5) 




-23.0 


70.9 (9.1) 


25.8 


54,725.82 


0.097 


-109.8 (8.8) 




-5.5 


124.1 (12.3) 


-0.1 


54,748.72 


0.347 


-123.8 (7.8) 




18.9 


189.0 (11.0) 


11.1 


54,754.67 


0.710 


106.2 (8.3) 




-17.5 


-216.7 (11.7) 


-21.6 
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Table 2 — Continued 



Date (HJD-2,400,000) 



(km s-i) 



Oi -Ci 
(km s-l) 



(km s-l) 



O2 -C2 
(km s-l) 



54,755.67 0.940 

54.757.65 0.395 

55.046.79 0.676 

55,070.71 0.161 

55.071.67 0.381 

55,072.71 0.619 

55.076.76 0.546 

55.078.66 0.982 

55,079.71 0.223 

54,341.82 0.752 

54.342.81 0.097 

54.343.82 0.449 

54.345.83 0.153 

54.347.82 0.847 

54.348.77 0.179 

54.397.65 0.250 

54,401.59 0.625 

54.402.66 0.999 

54.403.68 0.355 

54,696.81 0.723 

54.696.94 0.770 

54.697.80 0.069 

54.697.95 0.120 

54.698.74 0.397 

54,698.94 0.467 

54.699.83 0.777 

54.700.85 0.135 

54.701.86 0.488 

54,724.65 0.447 

54.724.75 0.481 

54.725.67 0.802 

54,725.77 0.837 

54.726.67 0.151 

54,726.77 0.185 

54.729.68 0.201 

54,729.74 0.226 

55.040.87 0.878 

55,041.67 0.155 

55.043.67 0.856 

55,044.77 0.240 

55,045.70 0.565 

55.046.68 0.906 



-61.1 (8.9) 


-33.1 


29 


3 


(12.5) 


11.9 


104.6 (10.7) 


1.7 


138 


7 


(15.1) 


11.7 


117.7 (12.9) 


-31.7 


-242 


5 


(18.2) 


-11.4 


133.8 (11.6) 


-5.1 


178 


7 


(16.3) 


20.2 


-97.8 (11.9) 


22.5 


143 


1 


(16.7) 


-3.5 


191.5 (10.5) 


7.4 


-232 





(14.8) 


47.7 


174.9 (13.8) 


4.2 


-249 


5 


(19.5) 


11.5 


-66.0 (10.3) 


-15.5 


109 


5 


(14.4) 


60.6 


141.3 (12.4) 


5.5 


219 


8 


(17.4) 


36.0 



MT771 



-139.3 (8.1 
-65.6 (4.0 
93.3 (6.9 
6.6 (6.3 
-154.9 (8.8 
23.7 (5.2 

70.3 (6.5 
3.3 (5.5 

-146.3 (4.9 

96.7 (4.1 
-79.5 (4.4 

-114.9 (5.7 
-81.1 (3.7 
-41.3 (5.5 
108.6 (5.8 
90.6 (5.6 
-125.0 (4.6 
-8.6 (26.3 
105.8 (5.4 
122.3 (3.8 
103.5 (4.4 
-138.8 (3.8 
-141.6 (4.3 
18.0 (3.9 
35.9 (3.1 
67.0 (3.6 

86.8 (2.7 
-123.0 (6.6 

11.9 (28.4 
-118.9 (6.8 

79.4 (9.0 
55.6 (6.4 

-123.4 (6.1 



-38.9 
-22.7 
-15.4 
-2.7 
-5.0 
-8.4 
-13.4 
2.9 
-23.2 
-23.7 
-0.4 
-2.7 
-12.9 
-19.6 
-11.4 
-11.2 
-8.3 
-1.8 
13.0 
13.0 
7.6 
-7.1 
4.9 
10.1 
-1.5 
17.0 
18.3 
33.2 
0.4 
33.4 
1.7 
9.0 
33.5 



85.1 (9.9 
19.3 (4.9 
-181.0 (8.5 
-69.2 (7.7 
113.7 (10.8 
-86.4 (6.4 
-144.7 (8.1 
-67.6 (6.8 
86.3 (6.0 
-179.5 (5.0 
91.6 (5.4 
109.7 (7.0 
56.0 (4.5 
26.9 (6.8 
-168.3 (7.1 
-147.3 (5.6 
106.6 (5.7; 

-125.6 (6.7; 
-146.6 (7.9 
-137.3 (5.5 
120.9 (4.7 

136.6 (5.2 
-60.8 (4.9 
-83.7 (3.8 
-70.7 (4.5 
-89.5 (3.3 
179.4 (8.2 

191.7 (8.3 
-107.9 (11.1 

-117.8 (7.9 
190.4 (7.5 
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Table 2 — Continued 





Vrl 


Oi 


-Ci 


Vr2 


O2 


-C2 


Date {HJD-2,400,000) 


(f) (km s^^) 


(km 




(km s^-"-) 


(km 


s-i) 


55,047.66 


0.250 103.3 (4.4) 




19.1 


-90.3 (8.9) 




36.8 



Table 3. Orbital Elements 



Element 


MT311 


MT429 


MT605 


MT696 


MT720 


MT771 


P (Days) 


5.752±0.005 


2.9786±0.0001 


12.27±0.01 


1.4692±0.0005 


4.3622±0.0003 


2.8635±0.0002 


e 


0.0 (fixed) 


0.0 (fixed) 


0.24±0.07 


0.01±0.02 


0.35±0.02 


0.05±0.03 


uj (deg) 


341±180 


339±180 


56±13 


170±180 


297±3 


41±29 


7 (km s^-*-) 


1±4 


-14±1 


-12±2 


-10±3 


-9±3 


-13±2 


To (HJD-2,400,000) 


52,149 (fixed) 


51,784.66 (fixed) 


54,653.2it0.5 


54729.0it0.7 


54,585.80±0.04 


54,408.4±0.2 


Ki (km s-l) 


88±5 


~140 


44±3 


261±4 


173±5 


139±5 


K2 (km s^l) 


118±9 


~205 


77±5 


307±9 


242±7 


164±6 


Ml (Mo) 


>3.0±0.7* 


7.6 


>1.3±0.3* 


>15.1±0.7'' 


>15.5±0.8a 


>4.5±0.3* 


M2 (M0) 


>2.2±0.4'' 


5.2 


>0.7±0.1'' 


>12.8±0.5^ 


>11.1±0.6'' 


>3.8±0.2^ 


/(m)i (Mo) 


0.4±0.1 




0.10d=0.03 


2.7±0.2 


1.9±0.2 


0.80±0.08 


/(m)2 (Mo) 


l.OiO.3 




0.5±0.1 


5.2±0.4 


5.3±0.6 


1.3±0.1 


S. Ci 


B2V 


[B3V] 


B0.5V 


09.5V 


BO-IV 


07V 


S. C.2 


B3V 


[B6V] 


[B2.5:V] 


B1:V 


B1-2V 


09V 


aisin i (R©) 


10.0±0.6 


8.3 


10.4±0.8 


7.6±0.1 


14.0±0.4 


7.867±0.004 


a2sin i (Rq) 


13.40±0.01 


12.1 


18.07±0.02 


8.921±0.004 


19.57±0.01 


9.3±0.3 


rmsi (km s^-'-) 


12.5 




7.2 


21.6 


17.9 


16.8 


rms2 (km s""*^) 


23.0 




8.6 


29.5 


25.4 


23.9 


in 


~42l' 


89 


~80'' 


~26'' 


~70'' 


<37'' 


Ti (K) 




19,100 










r2 (K) 




13,700 










Ri (R©) 




4.8 










(R©) 




3.6 










Note. — Calculated errors are located in parentheses. Square brackets indicate spectral types that have been determined indirectly 
from the data (i.e., based on the computed mass ratio rather than the component's spectrum). 

''Calculated mass is equal to Msin^i. 

''Estimated based on the theoretical masses of the components' spectral types. 
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Table 4. Ephcmcris for MT429 Spectroscopy 



Date (HJD-2,400,000) 





Vr 

(km s-i) 


51,805.89 


0.187 


0.8 (5.2) 


52,146.85 


0.658 


-28.0 (4.0) 


52,161.67 


0.634 


-23.1 (4.5) 


52,162.84 


0.029 


-5.1 (3.9) 


53,339.59 


0.100 


3.9 (5.7) 


53,340.59 


0.437 


-15.5 (4.5) 


53,989.78 


0.391 


-6.7 (6.6) 


54.286.91 


0.148 


2.2 (7.4) 


54,628.88 


0.957 


-15.9 (4.5) 


54,630.85 


0.619 


-22.8 (6.2) 


54,699.77 


0.758 


-32.0 (3.6) 


54,700.79 


0.099 


-3.7 (12.8) 


54,701.75 


0.423 


-12.7 (2.5) 


55,430.90 


0.221 


4.0 (7.6) 


55,497.59 


0.612 


-32.1 (6.5) 



Note. — The radial velocities listed in the table 
reflect the radial motion of the B2V component as 

diluted by the (apparently constant velocity) spec- 
trum of the brighter BOV component. 
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Table 5. OB Binaries in Cyg 0B2 



Star 


Type 


S.C. 


(days) 


c 


q 


Ref. 




SBl 


OSV -1- TRl 

vv o V 1 J 


4 8527 Cn 0002') 


n 1 1 en 04'! 


22-0 67^ 


X 


MT145 


SBl 


09III + [mid B] 


25.140 (0.008) 


0.291 (0.009) 


0.26-0.63=^ 


2 


MT252 


SB2 


B2III + BIV 


18-19 




0.8 (0.2) 


1 


MT258 


SBl 


OSV + [B] 


14.660 (0.002) 


0.03 (0.05) 


0.18-0.89=^ 


1 


MT311 


SB2 


B2V + B3V 


5.752 (0.005) 


0.0 (fixed) 


0.75 (0.07) 


3 


MT372 


SBl/EA: 


BOV + [B2:V] 


2.228 (fixed) 


0.0 (fixed) 


-0.6 


2,4 


MT421 


SB1:/EA 


09V + [B9V-A0V] 


4.161 




~0. 16-0.19 


5 


MT429 


SBl/EA 


BOV + [B3V] + [B6V] 


2.9786 (0.0001) 


0.0 (fixed) 


-0.68 


3,5,6,7 


MT605 


SB2 


BIV + [B2.5:V] 


12.27 (0.01) 


0.24 (0.07) 


0.57 (0.05) 


3 


MT696 


SB2/EW/KE 


09.5V + BOV 


1.4692±0.0005 


0.01 (0.02) 


0.85 (0.03) 


3,8 


MT720 


SB2 


BO-IV + B1-2V 


4.3622 (0.0003) 


0.35 (0.02) 


0.71 (0.06) 


1,3 


MT771 


SB2 


07V + 09V 


2.8635 (0.0002) 


0.05 (0.03) 


0.85 (0.05) 


1,3 


Schulte 3 


SB2/EA: 


06IV: + 09III 


4.7464 (0.0002) 


0.070 (0.009) 


0.44 (0.08) 


1,9 


Schulte 5 


SB2/EB 


07Ianfp + Ofpe/WN9 


6.6 (fixed) 


0.0 (fixed) 


0.28 (0.02) 


10,11,12,13, 






(+ BOV:) 








14,15,16 


Schulte 8a 


SB2 


05.51 + 06: 


21.908 (fixed) 


0.24 (0.04) 


0.86 (0.04) 


17,18 


Schulte 9 


SB2 


05: + 06-7: 


2.355 yr 


0.708 (0.027) 


0.9 (0.1) 


19,20 


Schulte 73 


SB2 


08III + 08III 


17.28 (0.03) 


0.169 (0.009) 


0.99 (0.02) 


2 


A36 


SB2/EA 


BOIb + BOIII 


4.674 (0.004) 


0.10 (0.01) 


0.70 (0.06) 


2,21,22 


A45 


SB2 


B0.5V + B2-3V: 


2.884 (0.001) 


0.273 (0.002) 


0.46 (0.02) 


2,22 


B17 


SB2/EB: 


07: + 09: 


4.0217 (0.0004) 


(fixed) 


0.75 (fixed) 


23,24 



Note. — Photometric types EW/KE, EA, and EB stand for eclipsing system of the W UMa type (ellipsoidal; P < 1 
day), Algol type (near spherical), and j3 Lyr type (e l lipsoi dal; P > 1 day) respectively. T he mass ratio for MT421 is 
calculated using the O star masses of iMartins et al.l | |2005| '| and interpolated AB masses of iDriUing &: Landolti ||2000|'I . 
The mass ratio associated with the system composed of MT429b + MT429c has been indirectly determined from the 
photometric data using PHOEBE. Square brackets indicate spectral types that have been determined indirectly from the 
data (i.e., based on the computed mass ratio rather than the component's spectrum). 

^The range in mass ratio incorporates both observational and inclination uncertainties 
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Fig. 1. — A composite of four spectra between a phase of = 0.457 and = 0.655 for MT311. 
The three panels show He I A5876 A, Ha, and He I A6678 A at quadrature respectively. 
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Fig. 2. — Spectrum of MT311 (top) obtained on 2001 August 24 and a composite of 
two spectra that best match the spectral types of the primary and secondary from the 
Walborn fc Fitzpatrickl (1l990l ) digital atlas shifted to the appropriate radial velocities. 
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Fig. 3. — Heliocentric radial velocity curve and orbital solution for MT311 using 11 of the 
highest S/N spectra. The filled points correspond to the primary (B2V) and the open points 
correspond to the secondary (B3V). The triangle represents an observation taken with Hydra 
at WIYN and the squares represent observations taken with WIRO-Longslit at WIRO. 
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Fig. 4. — Spectrum of MT429a (second from top) obtained on 20 01 September 8 and an 
09V (HD46202), a BOV (HD36512), and a BIV (HD144470) from the lWalborn fc Fitzpatri^ 
(1l990l ) digital atlas for comparison. 
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Fig. 5. — Comb i ned I -band photometric data from iPigulski fc Kolaczkowski 
Henderson et al.l (120111 ) folded at a period of 2.9786 days. 



(119981) and 
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Fig. 6. — Upper panel: the I-band photometric data and best-fit P HOEBE hght curve solu- 
tion fo r MT429. Open circles are the I magnitudes obtained from IPigulski &: Kolaczkowski 
(119981 ) , and the crosses are the I magnitudes obtained from [Henderson et al.l ( 1201 ll ) Lower 
panel: theoretical velocity curve for the putative B3V (solid curve) and B6V (dashed curve) 
components under the triple system (B0V-I-B3V-I-B6V) scenario described in the text for 
MT429. Data points (filled circles for WIRO, triangles for WIYN, a diamond for Keck) 
show the apparent observed heliocentric radial velocity in the composite spectrum which 
is dominated by the (presumably constant velocity) BOV modulated by the ~140 km s~^ 
amplitude signal of the fainter B2V. 
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Fig. 7. — Spectrum of MT 605 near quadrature (t o p) ob tained on 2001 September 9 and 
HD36960 (B0.5V) from the IWalborn fc Fitzpatricld (jl990l ) digital atlas for comparison. 
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Fig. 8. — Heliocentric radial velocity curve and orbital solution for MT605 using 10 of the 
highest S/N spectra obtained at WIRO with WIRO-Longsht. 
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Fig. 9.— He I A5876 A (left), Ha A6563 A (middle), He I A6678 A (right) in velocity space 
and in order of phase for observations of MT605 taken at WIRO between 2008 June 25 and 
2009 August 4. 
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Fig. 10.— He I A5876 A (left), Ha A6563 A (middle), He I A6678 A (right) in velocity space 
and in order of phase for observations of MT696. 
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Fig. 11. — Heliocentric radial velocity curves and orbital solutions for the eclipsing binary 
(W Uma type), MT696 (09.5V + BOV) using 17 spectra obtained at WIRO with WIRO- 
Longslit. The circles and solid line represent the velocities and orbital solution obtained 
from Ha and the triangles and dotted line represent the velocities and orbital solution ob- 
tained from He I A5876 APilled and open symbols correspond to the primary and secondary 
respectively. 
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Fig. 12. — A composite of seven spectra between a phase of = 0.5 and = 0.7 for MT720. 
The three panels show He I A5876 A, Ha, and He I A6678 A at quadrature respectively. 
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Fig. 13. — Heliocentric radial velocity curve and orbital solution for MT720 using 32 of the 
highest S/N spectra obtained at WIRO with WIRO-Longsht. 
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Fig. 14. — Heliocentric radial velocity curve and orbital solution for MT771 using 34 of the 
highest S/N spectra obtained at WIRO with WIRO-Longsht. 
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Fig. 15. — He I A5876 A in velocity space and in order of phase for observations of MT771 
obtained at WIRO with WIRO-Longsht from 2007-2009. 
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Fig. 16. — Map of OB stars in Cyg 0B2. Circles indicate the location of the 20 known 
binaries within Cyg 0B2. We find no apparent clusterings or groups of massive stars, nor a 
tendency for them to located near the center of the Association. 



